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Abstract 

We consider the d-dimensional diffusion equation dt<f>(x, t) — A0(x, t) 
with random initial condition, and observe that, when appropriately 
scaled, <fi(0, t) is Gaussian and Markovian in the limit d — > 0. This leads 
via the Majumdar-Sire perturbation theory to a small d expansion for 
the persistence exponent 0(d). Wc find 9(d) = id-0.12065...d 3/2 + . . . . 
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We consider the d-dimensional linear diffusion equation 

^ = A#(x, t ) (1) 

with initial condition </>(x,0) = ^(x), where ip(x) is a zero mean Gaussian 
random field of covariance (ip(x.)ip(x')) = <5(x — x'). We are interested in the 
persistence of the stochsastic process 3><f(i) = </>(0,i), that is, in the t — ► oo 
limit of the probability Q(t) that in the interval < t' < t the process 3>rf(i) 
does not change sign. This question was first asked by Majumdar et al. ||] 
and Derrida et al. f2j. The asymptotic decay of Q(t) appears to be given 
by a power law 

Q(t) ~ t-W (2) 

with a persistence exponent 9(d) whose precise value has been the main 
focus of study. One particular application concerns the coarsening of a 
phase-ordering system with nonconserved order parameter ||. 

Persistence exponents of random processes, also of interest in many other 
contexts in physics (see a recent review due to Majumdar Q), have turned 
out to be very hard to calculate. A perturbative method in this field was 
designed by Majumdar and Sire ||] and simplified by Oerding et al. f| 
(see also Majumdar et al. 0). It consists of expanding around a Gaus- 
sian Markovian process, for which the persistence exponent is known. This 
method was applied to zero temperature [|| and critical point || coarsening 
in the Ising model. More recently Majumdar and Bray [|| showed how to 
expand the persistence exponent of a Gaussian process in powers of a pa- 
rameter 1—p interpretable as a "fugacity" weighting the sign changes of the 
random process. 

For the diffusion equation (|T|) no exact expressions for 6(d) are known. 
Approximate numerical values come from simulation data Q , from the " in- 
dependent interval approximation" (HA) |IL 0, and from work by Newman 
and Toroczkai ||] based on a hypothesis concerning the "sign time distribu- 
tion" (see also Dornic and Godreche [[H]] and Drouffe and Godreche flTI). 
In this note we show that the persistence exponent 8(d) allows for an e ex- 
pansion in dimension + e. 



Using the Green function of the diffusion equation and writing for 
the surface of the d-dimensional unit sphere we can solve (|l|) and find 

o 1 / 2 poo . 

^ )= (4^ I drr^e-^VCr) (3) 

where is the appropriately normalized integral of the initial field f/K x ) 
on a spherical shell of radius r around the origin, 

tf(r) = S, 1/2 r -^ d -V lim -!- / dx-0(x) (4) 

Ar-+0 Ar J r <x<r+Ar 
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Hence ^(r) is Gaussian of mean zero and covariance (fy(r)^(r')) = 6(r — r'). 
We shall henceforth write d = e to indicate that the RHS of (||) is continued 
analytically to noninteger dimensions and that we prepare for the limit of 
zero dimension. 

The remaining discussion is easiest in terms of the new time variable 
r = e log t and the rescaled process 

$ £ ( T ) = (87r) £ / 4 e T / 4 $ e (e T / 4£ ) (5) 

which has the advantage of being stationary. Let us denote by Q(t) the 
persistence probability of 3> £ (t). This quantity then decays as Q(r) ~ 
exp(— 9(e)r) with e9 = 9. We consider now the explicit expression for the 
correlator F e (r — r') = ($ e (r)$ e (r)), which reads 

F e (T-T>) =( cosh ) (6) 

We observe that in the limit e->0we have the exponential decay Fq(t—t') = 
exp(— A|r — t'\) with A = \. It follows that $o(i) is Markovian and that it 
has 0(0) = A (see e.g. Q), whence we conclude that 9(0) = 0. 

For e > the process $ £ (t) is non-Mar kovian and there is no known 
way to calculate 9(e) exactly. We therefore consider F e — Fq as a small 
perturbation. This allows us to apply the principal result of the Majumdar- 
Sire perturbation theory Cast in a form first explicitly exhibited by 
Oerding et al. || and applied to the problem at hand it states that to 
lowest order in F e — Fq one has 

L vr Jo (1 — e _2Ar ) 3 / 2 J 

Upon substituting in (|7|) the explicit expressions for A, F e , and Fq, and 
multiplying by e, we get 



e r 1 [°°a (cosh£)- £ / 2 -e- T / 4 
4 L 1 ~ 2^ 7 T (1 _ e -r/2)3/2 



(8) 



We extract from (|8|) the first two terms of the e expansion by setting r = eu, 
expanding the integrand for small e, and performing an integration by parts. 
The result is 

a . . el / £ n3/2 roo du 1 

0{e) ~ 4-^2) J 7STT^ + --- 

8^)- 1 /2(i_V2) C(i) e 3/2 + ... 



c 
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= - - 0.12065... e 3/2 + ... (9) 
4 

Truncating this series after the e 3//2 term we find 9(1) = 0.1293... , which 
should be compared to the Monte Carlo estimate 9(1) = 0.1207±0.0005 M, 
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to the IIA result (9(1) = 0.1203... Q, @, and to the value 0(1) = 0.1253... of 
Ref. |J. No numerical values for 9(1) are available from the 1—p expansion 
of Ref. []|] . Finally, retaining in (pi) the full integral - which is the expansion 
proposed by the original authors ||] - gives for e = 1 the somewhat improved 
result 9(1) = 0.1267... . 

We conclude that the e expansion @ is significant as one of the rare an- 
alytical results that exist today for the persistence exponent of the diffusion 
equation. To make it numerically competitive higher orders in e should be 
calculated. 
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